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SUMMARY  OF  RESEARCH  ACCOMPLISHMENTS 


I.  Introduction  and  Background 

The  vibrational  potential  V  plays  a  central  role  in  interpretation  of 
diatomic  molecule  data.  Experimental  results  are  analyzed  in  terms  of  such 
potentials  and  the  potentials  in  turn  are  used  to  predict  unknown  properties. 
In  principle,  V  can  be  determined  by  solving  a  fixed-nucleus  eigenproblem  and 
correcting  for  nuclear-electronic  interactions.  In  practice,  however,  except 
for  1-  and  2-electron  systems,  such  ab  initio  theoretical  potentials  are  of 
quite  limited  reliability;  in  general,  ab  initio  Vs  cannot  be  used  to  deter¬ 
mine  accurate  vibrational  eigenenergies  E^  or  related  spectral  properties. 

As  a  result  of  the  intractability  of  the  many-electron  Schrodinger  equa¬ 
tion,  a  number  of  methods  have  been  devised  for  characterizing  and  relating  V 
and  Ey.  These  include  use  of:  power  series  expansions  for  V  as  a  function  of 
nuclear  separation  R,  approximate  integrals  from  which  V  is  found  at  a  finite 
set  of  R  values  through  inversion  of  vibration- rotation  data  (the  Rydberg- 
Klein-Rees  or  RKR  method),  long-range  multipole  forms  of  V(R),  and  closed-form 
expressions  for  V(R)  containing  adjustable  parameters. 

Unfortunately  each  of  these  methods  has  severe  limitations.  Near-equili¬ 
brium  properties  are  well  described  through  use  of  V(R)  power  series.  None  of 
the  methods  are  adequate  for  description  of  near-equilibrium  and  also  inter¬ 
mediate  energy  and  near-dissociation  properties.  Thus,  it  appeared  at  the  out¬ 
set  of  the  project  that  additional  research  should  focus  on  the  intermediate 
to  high  energy  range  of  diatomic  molecule  bound  states.  Better  understanding 
and  more  useful  techniques  were  both  sought.  Any  impro’/ements  achieved  would 
have  a  direct  beneficial  impact  on  Air  Force  laser  and  nuclear  weapons  effects 
technology  areas. 


The  proposed  research  had  two  specific  objectives: 


(a)  to  develop  expressions  for  the  vibrational  eigenenergies 

as  a  function  of  vibrational  quantum  number  v  useful  in  the  interme¬ 
diate  to  high  V  regions  for  diatomic  molecule  bound  states,  and 

(b)  to  evaluate  various  possible  causes  of  RKR  failure  near 
dissociation  and  to  devise  such  improvements  as  will  make  the  RKR 
procedure  useful  at  high  v. 

Consistent  with  these  objectives,  the  proposed  work  had  two  parts.  A  set  of 
model  vibrational  potentials  V  and  a  set  of  accurate  potentials  for  actual 
molecular  states  were  to  be  used,  each  with  a  variety  of  long-range  V  forms. 

For  these  potentials  the  eigenvalues  were  to  be  found.  To  aid  in  deter¬ 

mining  EyVS.  V  functional  forms,  reduced  mass  y  was  to  be  varied.  In  the  se¬ 
cond  part  of  the  work  the  effects  of  numerical  techniques,  low  density  of  data 
points,  and  the  breakdown  of  semi-classical  theory  at  high  v  were  to  be  studied 
for  their  roles  in  the  failure  of  RKR  near  dissociation. 

B.  Additional  Objective 

During  the  course  of  the  project  an  additional  objective  was  adopted: 

(c)  to  assess  the  accuracy  and  utility  of  rational  fraction  re¬ 
presentation  of  V(R). 

The  assessment  was  to  cover  a  variety  of  prototype  diatomic  molecule  elec¬ 
tronic  states. 

III.  Progress  and  Accomplishments 
A.  Computer  Codes 

Work  initially  centered  on  testing  available  computer  programs  for  eigen¬ 
value  extraction,  function  fitting,  and  RKR  inversion.^  As  a  result,  we  used 
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a  code  designated  SLEIGN  for  precise  eigenvalue  determination,  a  Householder 


reflection  code  for  linear  function  fitting,  a  Davidon  code  (called  VARMET)  for 
non-linear  fitting,  and  a  modified  RKR  program.  Both  VARMET  and  the  RKR  routine 
are  improved  versions  of  programs  that  we  originally  obtained  from  Kirtland  Air 
Force  Base  Computational  Services. 

SLEIGN  was  applied  to  the  Bom-Oppenheimer  (BO)  V  for  the  ground  state  of 

H2  to  obtain  precise  E^  with  both  the  true  reduced  mass  and  the  fictitious 

mass  9Uq.  This  potential  is  an  accurate  prototype  V  for  an  actual  molecular 
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state  with  long-range  -C^/R  behavior.  We  also  determined  first-order  WBK 
energy  levels  Ey(WBKj)  with  and  with  9y^  for  the  same  V. 

For  the  ground  state  of  H2  we  employed  SLEIGN  to  obtain  accurate  Born-Op- 
penheimer  E^  with  reduced  masses  y^  and  25  y^.  The  latter  high  mass  yields  77 
bound  states.  E  (WBK,)  were  also  determined  for  this  state  with  y  and  25  y  . 
Because  Ey(WBKj)  can  be  considered  a  continuous  function  of  v,  it  was  possible 
for  us  to  determine  WBK  eigenenergies  at  a  dense  mesh  of  noninteger  v  in  the 
immediate  vicinity  of  the  dissociation  limit  D. 

B.  Representation  of  E^  as  a  Function  of  v  (Objective  a) 

Ttere  is  theoretical  basis  for  a  power-series  expansion  of  E^  in  (v  +  5) 
near  equilibrium,^  i.e.,  at  low  v,  and  for  a  (v^^  -  v)"*  form  near  dissociation,^ 
i.e.,  at  high  v.  Here  v^  is  the  value  of  the  vibrational  quantum  number  at  the 
dissocation  limit  D,  and  m  is  a  parameter  dependent  on  the  long-range  form  of  V. 
Hence,  we  sought  a  composite  functional  form  of  E^  which  satisfies  these  two 
requirements  and  which  interpolates  well  at  intermediate  v.  In  particular,  we 
proposed  use^  of  the  expression  *  D-(Vp  -  v)"*  [L/N],  where  [L/N]  is  a  rational 
fraction  in  (Vp  -  v) .  The  degrees  of  numerator  and  denominator  polynomials 
are  L  and  N,  respectively. 

A  thorough  study  was  carried  out^  on  the  use  of  the  proposed  composite  E^ 
expression  for  the  H2*  ground  state  BO  eigenenergies.  A  variety  of  [L/N]  were 


investigated.  With  L  +  N  =  4,  a  polynomial  (N  =  0)  yielded  one  of  the  poorest 
fits.  As  our  proposed  form  is  nonlinear,  the  possibility  exists  of  having 
several  minima  in  each  space  of  variable  parameters;  hence,  it  was  necessary  to 
carry  out  least-squares  searches  from  a  variety  of  starting  points.  Our  best 
overall  fit,  with  L  =  2  and  N  =  4,  had  an  RMS  error  of  only  0.006  cm  ^  in  first 
differences  AE(v  +  s).  Difficulties  that  might  be  encountered  with  use  of  the 
expression  were  considered  in  some  detail.  The  effect  of  neglecting  terms  of 
degree  less  than  L  or  N  in  either  numerator  or  denominator  was  investigated. 

The  proposed  functional  form  for  E^  was  next  applied  to  the  H2  ground  elec- 
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tronic  state.  With  the  exception  of  [4/0]  all  L  +  N  =  4  fits  had  an  rms  error 
in  calculated  BO  AE(v  +  5)  of  0.02  cm  \  attributable  to  the  error  in  the  theo¬ 
retical  V.  The  error  in  calculated  D  -  Ej^^  (v  =  14  is  the  highest  bound  state) 
could  be  made  negligible.  Mass  25  levels  were  used  to  check  applicability  of 
the  proposed  expression  to  a  large  number  of  eigenenergies;  several  fits  yielded 
all  seventy-six  AE  with  an  rms  error  of  only  0.04  cm”^.  The  proposed  composite 
form  was  also  applied  to  the  experimental  vibrational  energies.  For  these,  best 
fits  had  an  rms  error  of  ^,0.07  cm  \  larger  than  for  theoretical  energies  pre¬ 
sumably  due  to  experimental  errors. 

The  same  form  that  had  been  applied  to  the  ground  electronic  state  of  H2 
was  used  with  isotopic  substitution  and  mass-reduced  quantum  numbers  to  predict 
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energy  differences  AE(v  +  j)  for  D2  and  HD.  Direct  conversion  from  H2  fittings 
yielded  rms  errors  from  experimental  AE's  of  0.49  cm  ^  and  0.43  cm  ^  for  D2  and 
HD,  respectively.  Least-squares  fitting  to  D2  and  HD  AE's  lowered  the  rms  errors 
to  0.034  and  0.045  cm  \  significantly  smaller  than  the  best  fittings  to  experi¬ 
mental  H2  data.  The  formula-predicted  dissociation  energies  differed  from  re¬ 
commended  theoretical  nonadiabatic  values  by  0.03  cm~^  for  D2  but  by  1.5  cm~^ 
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The  expression  for  was  applied  to  the  B  state  of  Br2,  a  prototype  with 
-Cg/R^  long-range  behavior.  We  fitted  the  functional  form  with  a  variety  of 
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[L/N]  to  the  fifty-five  experimental  first  differences  AE  for  Br2.  Our  best 
fits  each  had  an  overall  rms  error  of  0.015  cm"^  in  the  AE.  Apparently  erro¬ 
neous  experimental  results  were  exposed  through  consistently  large  errors  in  the 
same  AE's  among  different  fits.  Twelve  of  the  56  were  excluded  as  data  to 
test  the  interpolation  and  extrapolation  power  of  the  proposed  formula.  The 
rms  error  in  the  interpolated  AE's  was  essentially  the  same  whether  these  AE's 
were  included  or  excluded  as  fitting  data;  the  rms  error  in  the  extrapolated 

AE  at  the  three  lowest  v  increased  significnatly,  however.  Our  extrapolation 

-1  79 

for  D-E^^  yields  1.00  cm  .  Using  the  best  Br2  fits  and  mass  reduced  quantum 

81  - 1 

numbers,  we  fitted  to  Br2  and  found  an  rms  error  of  <0.03  cm  in  the  AE. 

Polynomial  fits  to  the  ^^Br2  E^  and  AE  were  attempted. The  results  indicated 

that  our  composite  functional  form  was  clearly  superior  to  polynomials  (a)  near 

the  dissociation  limit,  and  (b)  in  extrapolation  from  above  to  near- equilibrium 

energies.  Point  (b) ,  however,  requires  careful  scrutiny. 

127 

We  also  fitted  our  E^  vs.  v  representation  to  I2  B  state  experimental 
values^^  of  AE(v  +  J)  for  4  £  v  £  82.  Our  formula  proved  to  be  useful  for 
(a)  accurate  representation  of  known  energy  differences,  (b)  highlighting  and 
correcting  inaccurate  levels,  and  (c)  interpolation  and  extrapolation  of  missing 
I2  vibrational  energies.  Following  completion  of  our  calculations,  more  recent 
and  accurate  experimental  I2  data  were  brought  to  our  attention.  We  did  not 
have  opportunity  to  use  these  data;  however,  their  high  level  of  accuracy  and 
completeness  should  make  them  very  useful  for  testing  our  proposed  E^  formula. 

C.  RKR  Failure  near  Dissociation  (Objective  b) 

We  studied  in  detail  the  errors  in  RKR  turning  point  differences  6R  at  the 
highest-  v:»;r,iCional  levels,  v  =  60  for  H2^  (9  p^),  and  v  =  76  for  H2  (25  y^). 
With  a  cubic  spline  representation  for  E^  the  errors  for  the  SLEIGN  levels 
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and  were  41.1  a.u.  (atomic  units,  i.e.,  Bohr  radii)  and  12.9  a.u.,  respective¬ 
ly.  Even  for  the  WBK  levels,  the  errors  were  31.9  a.u.  for  E^q(WBKj)  and  7.4  for 
E^^(WBKp,  respectively.  The  latter  values  are  especially  significant  because 
RKR  determination  of  6R  is  a  first-order  WBK  inversion;  the  extremely  large 
errors  cannot  be  blamed  on  breakdown  of  the  semi-classical  theory.  Further,  our 
numerical  techniques  yield  very  precise  E^,  the  high  masses  lead  to  a  high  den¬ 
sity  of  data  points,  and  we  found  no  evidence  of  integration  difficulties  or 
related  errors.  We  attributed  the  error  to  the  E^  form  used  in  RKR  integrals, 
and  switched  to  our  proposed  form  described  in  Section  B.  When  we  represented 

Ey  =  [2/5]  (Vjj  -  v)^  in  the  RKR  integral  for  H2^(9  Uq)  >  the  6R  error  for  E^gCWBKp 

12  3 

was  reduced  from  31.9  to  0.0  a.u.  Similarly  with  [2/5] (Vg  -  v)  for  H2(25 

the  6R  error  for  E^g(WBKj)  changed  from  7.4  to  0.0  a.u.  The  error  reduction 

is  not  so  dramatic  for  R  associated  with  SLEIGN  E^;  nevertheless,  the  6R  error 

was  reduced  by  at  least  two  orders  of  magnitude  at  the  highest  integer  v  when 

(Vg  -  v)™  was  used  as  a  factor  in  the  E^  representation. 

Not  all  RKR  error  near  D  appears  attributable  to  the  E^  representation. 

However,  incorporation  of  the  long-range  term  (Vg  -  v)"*  in  a  hybrid  representa- 

4  6 

tion  for  E^  removes  the  key  source  of  error  at  high  v  both  for  1/R  and  1/R 
prototype  states.  It  is  not  clear  that  significant  additional  effort  is  war¬ 
ranted  in  isolating  causes  of  residual  error. 

D.  Rational  Fraction  Representation  of  V  (Objective  c) 

In  work  not  originally  proposed,  we  fitted  rational  fractions  [L/N]  in  R 

to  V(R)  for  the  ground  states  of  H2*;  L  and  N  were  chosen  to  give  the  correct 

13  14 

long-range  form  for  V.  A  [5/9]  fitting  yielded  ’  an  approximate  V  with  a 

maximum  error  of  0.04  cm  ^  in  the  domain  1.0  <  R  <  100  a.u.  This  V  had 

—  —  approx 

only  a  single  pole  for  real  R  2^  0  and  this  was  at  0.002  a.u.  (This  pole  should 
be  exactly  at  R  =  0.)  Further,  all  E^  but  E^g  obtained  by  solving  a  Schrodinger 
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equation  containing  this  approximate  V  are  accurate  to  within  0.01  cm”^  and 
is  only  in  error  by  0.016  cm"^. 

The  success  of  the  rational  fraction  potential  representation  led  iis 'to 
apply  this  expression  to  the  Bom-Oppenheimer  potential  for  the  ground  elec¬ 
tronic  state  of  H2.  We  restricted  N  to  L  +  6  to  yield  the  correct  long-range 
form  for  this  state.  Our  best  [5/11]  fit  had  an  rms  error  from  Kolos  and 

Wolniewicz'  V  of  0.01  cm  ^  on  0.7  £  R  £  12  Bohr  radii. The  [5/11]  fitting  has 
-4 

a  pole  at  1.9  X  10  a.u.  with  a  residue  of  0.98  instead  of  the  correct  values 
of  0  and  1,  respectively. 

Next,  we  applied  the  V(R)  expression  to  the  2pTiy  state  of  This  state 
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is  bound  but  has  a  repulsive  long-range  1/R  form.  Thus,  with  a  potential  mini¬ 
mum  at  R  =  7.9  a.u.  and  a  small  maximum  at  R  "v  26  a.u.,  this  state  places  strong 
demands  on  a  V(R)  representation.  In  spite  of  this,  we  found  that  a  [6/9] 
fitting  achieves  an  rms  error  at  74  R  values  in  5  £  R  £  90  a.u.  of  only 
0.0014  cm  ^.^^The  [6/9]  is  reliable  for  interpolation.  It  extrapolates  at  low 
R  to  a  pole  at  -1.5  x  10  ^  a.u.  with  residue  0.92  in  place  of  correct  values 
0  and  1,  respectively.  The  reliability  of  the  extrapolation  near  R  =  0,  even 
though  no  fitting  points  below  5  a.u.  were  used,  is  probably  due  to  the  domi¬ 
nance  at  intermediate  R  of  the  nuclear  repulsion  over  the  purely  electronic 
energy. 

Rational  fraction  representation  applied  to  the  Bom-Oppenheimer  potential 
V  for  H2*  2pay  is  of  interest  because  this  is  a  prototype  -1/R^  van  der  Waals 
electronic  state.  A  [6/10]  fraction,  least-squares  fitted  on  0.8  <  R  <  90  Bohr 
radii,  has  an  rms  error  of  0.003  cm"^  at  fitted  points. The  [6/10]  accurately 
reproduces  the  shallow,  large-R  minimum,  and  yields  reasonably  reliably  the 
Coulomb  pole  at  R  =  0. 

The  vibrati'  pot*"  lal  of  the  3dJg  electronic  state  of  H2*  was  renresented 
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by  rational  fractions  in  nuclear  separation  R  over  the  range  R  =  1  to  90  a.u. 

2 

This  state  with  a  deep  well  is  characterized  by  a  strong  1/R  attraction  with 
significant  higher-order  raultipole  contribution  even  at  very  large  R.  The  long- 


range  form  of  the  potential  was  reproduced  by  using  either  one  or  three  con¬ 
straints  on  the  denominator  polynomial.  The  best  overall  fit  achieved  over  the 
entire  range  (1-90)  has  an  rms  of  0.03  cm  However,  fits  over  restricted 

ranges,  1-28  a.u.  and  5-90  a.u.  each  achieve  an  order  of  magnitude  smaller  error. 


To  complete  consideration  of  prototype  V,  three  purely  repulsive  states  of 

H2*  were  considered.  These  are:  the  ^dir^  with  1/R^  long-range  repulsion  and 

2 

promotion  from  n  =  2  at  R  =  «>  to  3d  at  R  =  0,  2sOg  with  1/R  repulsion  and  no 
p?.'omotion,  and  Spo^  with  1/R“  repulsion  and  promotion  from  n  =  2  to  3p  in  the 


united  atom.  In  spite  of  this  diversity,  all  three  had  best  fits  on  1  £  R  £  90 
a.u.  with  rms  errors  of  about  0.015  cm  Work  on  these  three  states  has  not 


yet  been  reported  at  a  scientific  meeting  or  submitted  for  publication  in  a 


journal. 

Rational  fraction  representations  have  thus  yielded  a  very  high  level  of 
accuracy  for  eight  diatomic  states  with  diverse  V(R)  forms.  Generally,  a  single 
fraction  was  applicable  over  a  very  wide  R  range. 
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